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ABSTRACT 

In the course of the formation of cosmological structm'es, large shock 
waves are generated in the intra-cluster medium. In analogy to processes in 
supernova remnants, these shock waves may generate a significant population 
of relativistic electrons which, in turn, produce observable synchrotron emis- 
sion. The extended radio relics found at the periphery of several clusters and 
possibly also a fraction of radio halo emission may have this origin. Here we 
derive an analytic expression for (i) the total radio power in the downstream 
region of a cosmological shock wave and (ii) the width of the radio-emitting 
region. These expressions predict a spectral slope close to -1 for strong shocks. 
Moderate shocks, such as those produced in mergers between clusters of galax- 
ies, lead to a somewhat steeper spectrum. Moreover, we predict an upper limit 
for the radio power of cosmological shocks. Comparing our results to the ra- 
dio relics in Abell 115, 2256, and 3667, we conclude that the magnetic field in 
these relics is typically at a level of 0.1 /nG. Magnetic fields in the intra-cluster 
medium are presumably generated by the shocks themselves, this allows us 
to calculate the radio emission as a function of the cluster temperature. The 
resulting emissions agree very well with the radio power-temperature relation 
found for cluster halos. Finally, we show that cosmic accretion shocks gener- 
ate less radio emission than merger shock waves. The latter may, however, be 
detected with upcoming radio telescopes. 

Key words: galaxies: clusters: general - radio continuum: general - diffuse 
radiation - radiation mechanisms: non-thermal - methods: analytical 



1 INTRODUCTION 



Diffuse radio objects that extend over several Mpc but have no optical counterpart have 
been observed in several clusters of galaxies. These objects fall into two categories: radio 
halos are centrally located in the cluster, they are fairly regular in shape and show very little 
polarisation. In contrast, radio relics appear at the periphery of clusters, have a filamentary 
morphology and are often highlv p olarised. Abell 3667, for instance, shows a spectacular 
double relic (|Rottgering et al.lll997l ). where the two radio objects, separated by ~ 4 Mpc, 
straddle symmetrically the cluster X-ray emiss ion. Simil arly extended and remo te relics have 
also been found in All 5 (" Govoni et all 120011 ). A2256 (|Giovannini et al.lll999h . and A2345 



( Giovannini et al'IE'999. ) to name just a few. All of these clusters also show sign s of an ongoin; 



or recent merger: A115 has a double peak in the centre in X-ray observations ()Shibata et al 



a 
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1991 . A2256 shows also a double peak (|Sun et al.ll2002[ ). A2345 is ass umed to be a dynami- 



cally young system since it shows multiple X-ray substructures (see discussion in lDahle et al. 
I2nn2| ). and in A3667, among other indications, a cold front h as been found which is presura - 
ably related to a slightly supersonically moving substructure fIVikhlinin fc Markevitchll2nn2h . 
In particular, the last example sugg ests that the re lics are located at the shock wave that has 
been produced by the substructure. iRoettiger et al. (1999) were able to reproduce the X-ray 
morphology and the location of the radio emission for this particular cluster in a merger 
simulation. 

Even if the X-ray emission at the periphery of galaxy clusters is too faint to identify 
clearly shock fronts, radio relics seem to trace cosmological shock waves. Two models have 
been proposed for the origin of radio relics: They could ei ther be caused by diffusive shock 



acceleration of electrons at the shock waves themselves ( Enfilin et al. 19981 : Miniati et al. 



20011 1 , or they could be old radio bubbles that are compresse d by the passing shock wave and 



thus induced to em it obser vable synchrotro n emission again ( Enfilin fc Gopal-Krishnal 2001 



lEnfilin fc Briiggenl [2002: Hoeft et al.ll2004l ). The former rehcs are sometimes called radio 
gischt, while the latter ones are known as radio phoenix ( Kempner et al. 2004[ ). The radio 
relic in A85 ( Slee et al may serve due to its size, morphology, and strong polarization 

as prototype for the radio phoenix class. 

Radio halos, in contrast, do not directly trace the shock front s. Instead, the radio mor- 
phology is similar to the overall X-ray morphology o f the cluster (iGovoni et al.ll2001 ). The 
radio spectrum of halos tends to steepen with radius ( Feretti et al. l2004l2005f l."This may be 
regarded as evi dence that the emitting electrons are accelerated by turbulence in the ICM 
()Brunettill20"ol . which is a result of a merger. However, the origin of the halos is still unclear. 
There is increasing evidence that halos are indeed related to cluster mergers, but it is a sub- 
ject of a debate whether the relativistic electron population that is responsible for the radio 
emission is emitted by a central AGN, whether it is produced by strong shocks, whether it 
stems from the decay of relativistic protons or whethe r it is genera ted by turbulence in the 
downstream region of the shocks (for a review see e. g. lFerettill200,^. Interestingly, the radio 
power of halos sc ales with the cluste r X-ray himinositv (jFerettill200,^ or, equivalently, with 
its temperature ( Keshet et al.ll20o3 ). Therefore, the non-thermal component in a galaxy 
cluster, if existent, seems to be related to the thermal plasma. 

Numerical simulations show that the formation of the large-scale structures leads to 
a variety of shocks in the intra-cluster medium (ICM) and inter-galactic m edium (IGM) 
()Miniati et al.ll2000l : IRvu et al.lE)oi IPfrommer et"^l200fil : iHoeft et ahlboOfih . These shocks 
are crucial for the thermal state of the ICM. As matter falls into the deep potent ial wells of 



galaxy clusters, its shocks to the high temperature of the ICM of ~ 10^ — 10^ K. iRvu et al 



f)2003[ ) distinguish between 'internal' and 'external' shocks. The former affect material that 
is already heated to roughly the cluster temperature, whereas the latter heat rather cold gas 
for the first time. Hence, internal shocks are typically weak while external shocks are very 
strong. 



Cosmic rays with energie s up to ~ lO^^eV are 



of supernova remnants (e.g. iBerezhko et al.ll200 



resumably accelerated at the shock fronts 
Vink fc Lamind 20031 ). The acceleration 



is believed to be caused by diffuse shock acceleration: Particles may be scattered repeatedly 
between the upstream and the downstream regions, separated by t he shock. Each tirne the 
particl e crosses the front, it gains ki netic energy (see for a review lAxford et al.1ll978l : Isdl 



Il978al f9: Blandford fc Ostriker 19781 ). Even though the physical conditions of shock waves 



in the ICM are somewhat different from those in supernova remnants, one would also expect 
diffuse shock acceleration to be efficient at cosmological shocks. 

Various groups have tried to estimate the radio emission of the shock fronts by combining 
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the probability distribution of cosmic shock fr onts with m odels for diffuse shock accelera- 
tion and the subsequent synchrotron emission. iKeshet et al. (2004) found that 10 % of the 
extragalactic radio background below 500 MHz should be caused by cosmic structure for- 
mation shocks for ACDM cosmology. They assumed that the emission comes mainly from 
external accretion shocks, in particular those that surround massive clusters of galaxies. 
Essential for the radio emission is the strength of the magnetic field since the emission is 
synchrotron emission. Unfortunately, the magnetic field distribution in clusters of galaxies 
and their surrouiidings is difficult to determine and still poorly constrained (for a review see 
Govoni fc Ferettill2004^ . 



In this paper, we compute the radio emission in the downstream region of a shock front. 
We assume that diffuse shock acceleration produces a relativistic electron population with a 
power-law distribution in the energy spectrum (Sec. l2.2j) . As the plasma moves downstream, 
the high-energy electrons cool by inverse Compton and synchrotron losses (Sec. l2.3j) . The 
main objective of this paper is to derive an analytic expression for (i) the total radio power in 
the downstream region of a cosmological shock wave and (ii) the width of the radio-emitting 
region. In particular, we wish to explore the dependence of the power, spectral shape and 
spatial extent of the radio emission on the magnetic field and the shock parameters. 

We apply our model to the relics observed in Abell 115, 2256, and 3667 and infer the 
magnetic field strength in the relic region (Sec. IS}- Assuming that the magnetic field energy 
density is a fixed fraction of the thermal energy allows us to express the radio power as 
a function of the cluster temperature (Sec.|3I). Finally, we make predictions for the radio 
emission of cosmic accretion shocks (SecE}. 



2 RADIO EMISSION BY NON-RADIATIVE SHOCKS 
2.1 Non-radiative shocks 

A shock surface separates two regions: The upstream plasma moves with velocity towards 
the shock front, while the downstream plasma departs with Vd- As the plasma passes through 
the shock front, mass, momentum, and energy fiuxes are conserved, which is expressed in 
the Rankine-Hugoniot relations: 

PuVn = PdVd 

Pu + PuVl = Pd + Pdvl (1) 

1 2 Pu 1 2 Pd 

2 Pu 2 Pd 

where p denotes mass density, P pressure, and u specific internal energy. On timescales 
relevant for the shock propagation, we assume that the intra-cluster medium behaves like a 
polytropic, ideal gas. Hence, the pressure can be written as 

P=(7-l)pM, (2) 

where 7 is the adiabatic index, and the specific internal energy depends only on the gas 
temperature, T, 

(7 - l)pmp u = ksT, (3) 

where rup is the proton mass and p is the molecular weight. We determine the latter by 
adopting a fully ionised plasma with primordial chemical composition. In contrast to the 
conserved properties, the entropy of the plasma is increased by the dissipation at the shock 
front. For simplicity the term 'entropy' will refer throughout this paper to a monotonic 
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function of it, namely the entropic index, 

S = up^-\ (4) 

The strength of non-radiative shocks in a polytropic gas can be characterised by a single 
parameter, e. g. the compression ratio 

r = -, (5) 
or, equivalently, by the entropy ratio 

d - ^- (6) 

The ratio of the specific internal energies can be given as a function of the compression and 
entropy ratios 

— = ^ ( — ] =(ir^ ■ (7) 

Combining the conservation laws, Eq. (jSj), and using the definition of the entropy, Eq. 
allows us to relate the compression and the entropy ratios by the implicit equation 

(7-l)grT+(7 + l) ^ ^ 

Hence, if we can determine the entropy ratio, q, for a shock we can compute the compression 
ratio by Eq. (jH| and finally the ratio of internal energies by Eq. ((7)). 

Customarily, the strength of a shock is characterised by the upstream Mach number 

-M = (9) 

where Cu denotes the upstream sound speed, which depends on the specific internal energy 

by 

c2=7(7-l)n„. (10) 
Invoking flux conservation, Eq. (0), we can rewrite the Mach number 



M' = -P^ = - (11) 

cl Pn Pd- Pn 7 r -1 



For large Mach numbers the entropy ratio is proportional to 7W^, see Fig.^ For later use 
we rewrite the downstream velocity 

2 7-1 qr^ -1 _ ^2 f-,o\ 

Vd = ; r = Md, (12) 

qr^ r — I 



2.2 Diffuse shock acceleration 



In supernova remnants, there is evidence that electrons and protons are accelerated by diffuse 
shock acceleration (DSA) to energies of ~ 10^^ eV (e.g. Berezhko et aLlbOO^I : Vink fc LamingI 
2003| ). In DSA, particles are accelerated by multiple shock crossings, in a first-order Fermi 
process. If the shock thickness is much smaller than the diffusion scale which in turn has to 
be much smaller than th e curvature of the shock front a o n e- dimensional diffusion-conv ection 
equation can be solved (lAxford et al.lll978l : lBdlll978al lbl: iBkndford fc Ostrikeilll978h . The 
result is that the energy spectrum of suprathermal electrons is a power-law distribution, he (x 
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E ^. The spectral index, 
ratio at the shock front 
r + 2 



s, of the accelerated particles is only related to the compression 



(13) 



r — 1 

For strong, non-radiative shocks with Mach number ^10, the compression ratio is always 
close to 4, hence the slope s is always close to 2. This allows to explain the spectrum of 
cosmic rays over a huge range of energies and may b e cons idered as a piece of evidence 
for PSA. For a review of diffu s e shock acceleration se e Drurv (198i|); Blandford &: Eichlei 
( 19871 ): IJones fc EUisonl (|l99lh : iMalkov fc O'C DruT^ ^20Ql\) . 

Diffuse shock acceleration theory suffers from the complexity of real shock fronts and 
the highly non-linear interaction between different processes such as cosmic rays and mag- 
netic waves. The efficiency of accelerating electro ns and protons has been inferred from the 
observation of supernova remnants. For instance, iDver et al.l (|2001l l estimated for SN 1006 
that a few percent of the shock energy is transferred to supra-thermal particles. Similar to 



Keshet et al. 



( 2004[ l. we assume that a fixed fraction of thermal energy injected at the 
shock front goes into the acceleration of suprathermal electrons. 

Moreover, we multiply the spectrum by an upper cut-off factor since electrons can only 
be accelerated to finite energies. Typically, the maximu na energy is estimat ed by comparing 



the acceleration e-folding time and the cooling time ( Keshet et al. 2003| ). We use here a 



smoothed high-energy cut-off, the reason for this particular choice will be become clear in 
the next section. Thus, the electron spectrum generated by DSA, is given by 



UEiE) 



dE 



a 



spec 



s-2 







6 < Cmax 

elsewhere 



(14) 



where we have used the abbreviation e = E/rUeC^. The normalisation constant can be 
interpreted as follows: We compute the electron number density of electrons with energies 
in the interval [E, E + AE]. If we choose E = AE = rrieC^, the number density becomes 



UeAE = UeCi 



spec- 



nipC 



Hence, Cgpec gives basically the fraction of electrons at E 

We determine the fraction of suprathermal electrons by postulating that a fixed frac- 
tion of the energy injected at the shock front goes into the acceleration of electrons to 
suprathermal energies. After passing the shock, the downstream gas has gained a thermal 
energy of 

MdPd - r^u^Pu = UdPd , 

Q 

where r^u^p^ would be the energy density after pure adiabatic compression. Hence, the 
spectrum can be normalised by 

1 



dE ue^E) E = ieUdPd 



(15) 



A crucial parameter is the minimum energy, E^[^, above which electrons are considered to be 
suprathermal. It is particularly important for spectra significantly steeper than s = 2 since 
for these spectra most of the suprathermal energy is carried by electrons immediately above 
-Emin- Therefore, the normalisation of the spectrum depends on ii^min- Hybrid simulation of 
a collisionless shock indic ate that there is continu ous transition from the thermal to the 
suprathermal distribution ( Bennett fc EUisonl 19951 ) . Simulations of electron acceleration at 
high Mach number shocks may rely on injection of electrons with a fixed energy. However, the 
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typical injection energy corresponds to thermal energy in the ICM ( Levinsonlll994[ l . therefore, 



the suprathermal spectrum is expected to be a continuous extension of the thermal spectrum. 
We assume here that the thermal Maxwell-Boltzmann distribution goes over continuously 
into the power-law spectrum of the suprathermal electrons 

^Maxwell ^^^^^^^ = ^^(Einin). 

This condition leads to an implicit equation for -Emin that has to be solved simultaneously 
with the normalisation of the spectrum. Basically, the transition energy is tightly coupled 
to the temperature of the plasma, 

E^in ~ 10 X ^bT, 

if no extreme parameters are chosen. 

We can write for the normalisation of the electron spectrum 

q-l 2 ~ ~l-s \ ^ 2 

^dPd = C*spec UleC I dc 6 * < 1 — > = He C*spec f^eC Ispeci 

where 



oo 



s-2 



/spec = / de e'-' <^ 1 - ^ ^ . (16) 

J Bmin L 6 max J 

Hence, the constant becomes 
^ _^Udmp (g- 1) 1 

L'spec Se 2 r ' \ / 

^ C ^ TTLq^ (J 1 spec 

Cfpoc Cipcc 

where the first factor, C^pec, includes the plasma contributions, while the second factor, C^^^^, 
incorporates the dependencies from the shock strength. Since the electron spectrum is steep 
for Mach numbers < 2, the fraction of relativistic electrons is low in this regime, see Fig.|21 
Note that also the temperature of the downstream plasma affects the relativistic electron 
density. For Mach numbers > 10 the fraction of relativistic electrons is virtually constant 
since the electron spectrum varies only negligibly. 

In this paper, we only consider the acceleration of electrons in order to compute the 
radio emission from cosmological shocks. Clearly, there is evidence that also protons are 
accelerated to highly re l ativis tic energies at shocks, probably even more efficiently than 
electrons, c.f. Levinsoiil ( 1994 ). Because of their larger mass, protons cause significantly 



less synchrotron emission than electrons. However, in inelastic collisions, relativistic protons 
produce electron/positron pairs which, in turn, radiate synchrotron emission. In the ICM, 
though, the life time of relativistic protons is about the Hubble time ( Denniso'nl 198Cl( ). As a 



result, relativistic protons gradually accumulate in the cluster potential. While the emission 
by secondary electrons may contaminate to some degree the emission by primary electrons, 
it is unlikely to account for the bulk of the radio emission from radio relics. Radio relics 
have a filamentary morphology and the emission is confined to regions close to the shock. 
This morphology is difficult to explain with superthermal protons and is more likely caused 
by short-lived, primary electrons. The emission from secondary electrons and positrons pro- 
duced by relativistic protons may be responsible for the cluster-wide emission observed 
in radio halos. However, acceleration of protons has another consequence. Suprathermal 
protons can affect the pressure in the gas and thus affect the structure of the shock. An 
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admixture of relativistic protons lowers the adiabatic exponent, 7, to below the value for a 
non-relativistic gas of 5/3. In the extreme case of a fully relativistic gas, 7 is 4/3. An effect of 
a lower 7 is that the compression ratio of the shock can be higher, and, as a result, Eq. ()13|1 
becomes inaccurate. Moreover, the back-reaction of the relativistic particles on the shock 
leads to non-linear particle acceleration, which results in flatter particle spe ctra than the 
power-laws obtained in the test-particle approximation. Amato fc Blasil ( 2005t ) have demon- 
strated the effects of non-linear particle acceleration on shock structure and particle spectra. 
While there is some evidence for non-linear effects in supernov a remnants, the details of this 
process remain controversial ( Ellison &: Cassam-Chenai2005^ . In particular, it is unclear to 
what extent non-linear shock acceleration takes place in cosmological shocks. Here we will 
adopt the results obtained with the linear test-particle approximation. 



2.3 Radiative cooling of high energy electrons 

High-energy electrons in the ICM cool mainly by inverse Compton scattering (IC) with CMB 
photons and by synchrotron radiation losses 

^ = (18) 

at rUeC^ 

where E is the energy of an electron. The cooling constant is given by 

4 (Tx 

Ccool = 7; {mcmb + Ub}, (19) 

where ctt is the Thomson cross section, mcmb and ub are the energy densities of the radiation 
background and the magnetic field, respectively. We have neglected the effects of adiabatic 
cooling and heating which depend on the local flow conditions. Adiabatic cooling can in 
principle be included but for our purpose we will ignore it here. Electron number conservation 
leads to a PDE for the energy spectrum 

-g^ = Ccooie + 2Ccooie nE, (20) 
with the solution ( Kardashev 1962t ) 



nE{E,t) = { ^^'^'^P'^^^i^^ '{^~ + • e Ccool t < 1 - e/e„iax (21) 

: elsewhere 

We have chosen the cut-off factor in Eq. p4|) consistently with this solution, hence the initial 
spectrum is already a solution to Eq. (j^ . 

The energy densities of the cosmic radiation background with temperature Tcmb = 
TcMB,o(-2 + 1) and of the downstream magnetic field, 5d are 

WCMB = «rad^CMBO (^ + '^Y = ~S (2^) 

= (23) 

respectively, where z is the redshift and Orad the radiation constant. For convenience, we have 
introduced a magnetic field equivalent for the background radiation, -Bcmb- For z = 0, Bqmb 
is about 3.24 /xG. Since we are primarily interested in radio emission which is produced by 
relativistic electrons, it is convenient to express the related equations in terms of the Lorentz 
factors, 7l, of the electrons. Only electrons with 7l ^ 1 contribute significantly to the radio 
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emission and we can approximate the Lorentz factor 
E E 

7l = ^ + 1 ^ ^ = e for E > m^c^ (24) 

For any observing frequency, there is a minimum Lorentz factor below which the radio 
emission is negligible. On the other hand, Eq. ()21|) shows that the electron spectrum has a 
time-dependent, maximum energy. For any Lorentz factor we can determine the time when 
then maximum energy in the spectrum equals the given 7l, which is given by 

tcool(7L) = ^ \i . = 2.32 X 10^2 + i)-4y^^ 

where the right expression is computed for pure inverse- Compton cooling (we have assumed 
here e/cmax 1). After a time tcooi(7L) there are no electrons with Lorentz factor above 7l. 



2.4 The synchrotron emission 

The synchrotron power of a single electron , Pp, that moves with a L orentz factor 7l in a 
magnetic field B is (for an introduction see Rvbicki fc LightmanI 19861 ) 



d-Pe(7L, t^obs) _ VSe^Psina ^ / z^bsA ^ ^ ^ f 

dv m^c^ \ i^c J ~ V / ' 

where e is the electron charge, a the angle enclosed between the direction of the motion and 
the magnetic field, and Uc is the characteristic frequency 

S'-yfeB sin a 
An rUeC 

In the following, we assume that the directions of the electron motion and the magnetic field 
are uncorrelated, therefore we take the average value (sin a) = 7r/4, where the average is 
over solid angle. The function F{x) has to be integrated numerically. 



Fix)=x / d^K.iO 

J X 



where K is the modified Bessel function. The emission per volume is given by the convolution 

12 T3(,, \ /"oo /"oo 



d^P(^^obs) 



dVdu 

We introduce a new variable 



POO POO 

/ dE He Pe = Csync / d7L rUeC^ UEilL, t) F(7l, Z/obs) 

Jo Jo 



1 1 

2 



..y^^,,.au^4.85xlO-|4)'|H^y ... (26) 

V i^obs 16 rrieC [l^^ ) I t'obs J 



Substituting 7l by r, the integration becomes 

d^P(z/obs) _ CsyncirieC^ 



dv du a 



^ dTnEiT,t)Fi-\. 



(27) 



Assuming a power-law spectrum for he, we compute the cumulative emission, see Fig.Ol 
The bulk of the emission comes from Thaif ~ 0.03. This implies that electrons with Lorentz 
factors about 



1 ^ 



-,.«^6.2X10M^ l^j' (28) 
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contribute most to the synchrotron emission. Note that at 7l ~ 3 x 10^ and below, Coulomb 
losses coristitu te the most important cooling mechanism for suprathermal electrons, see e.g. 
barazi Hence, at very low observing frequencies in regions with a rather strong 



magnetic field. Coulomb cooling becomes important. For instance, Lofar may observe at 
30 MHz a radio halo region where the magnetic field is of the order of 3/iG. Eq. ()28p indicates 
that the bulk of the radio emission would come from electrons with Lorentz factor 7l ~ 
50. For Lorentz factor that low Coulomb cooling is important, hence, our analysis may 
overestimate the synchrotron emission at very low observing frequencies. 



2.5 Total emission behind the shock front 

In Sec. 12. 21 we worked out the relativistic spectrum of electrons accelerated at the shock front, 
in Sec. 12. 31 we worked out the evolution of the spectrum subject to cooling, and in Sec. 12. 41 
we computed the synchrotron emissivity for a homogeneous electron population. Now we 
combine these results, to compute the total emission behind the shock front. We adopt 
the following scenario: Diffuse shock acceleration generates a population of suprathermal 
electrons. In energy space, these electrons are distributed initially according to a power-law. 
While they advect with the downstream plasma, the relativistic electrons cool by synchrotron 
and inverse Compton losses. With increasing distance to the shock front, the maximum 
electron energy decreases and the radio emission diminishes. The total emission is obtained 
by summing up all contributions from the plasma from the shock front to the distance where 
the electron spectrum is too cool to allow any further radio emission. 

To compute the spatial distribution of the emission, we have to model the downstream 
region of the shock front. Instabilities of waves propagating at the shock surface presumably 
originate turbulence in the downstream plasma immediately behind the shock surface. If this 
turbulence can also be found at larger distances from the shock s urface depends stro ngly on 
the viscosity of the ICM, which is still rather uncertain (see e. g. iFabian et al1l200,5h . If the 



viscosity is indeed close to the Spitzer value of a hot plasma, as suggested by the Ha fila- 
ments surrounding NGC 1275, the formation of turbulence would be significantly hampered. 
Hence, we assume the flow in the downstream region to be steady and laminar. 

Diffusion can alter the spatial distribution of suprathermal electrons. In the absence of 
magnetic fields, the diffusion length is dominated by Coulomb collisions. The presence of even 
very weak magnetic fields suppresses diffusion in various ways. In directions perpendicular to 
magnetic field lines electrons propagate roughly with the Bohm diffusion rate. The resulting 
distance covered within the cooling time of even very energetic particles is negligible for our 
purposes. Parallel to the magnetic field lines the electrons propagate with the drift velocity, 
provided the gyroradius is significantly smaller than the coherence length of the magnetic 
field. The gyroradius of relativistic electrons is given by 

mpC^ „ « , f 1l\ f B 



7l — — ?a 6 X 10 kpc 



eB ' VlOV yifiG^ 

Hence, the gyroradius of suprathermal electrons in the ICM is of the order of 10^^ — 
10^^ cm. The correlation length of magnetic fiel ds generated at the shock surface is expected 



to be of the order of 10 cm (jMedvedev et al.. .2006>) . but it is also expected to increase 



rapidly ( Medvedev et al. 2005t ). As a result, electrons drift virtually freely along the field 



lines. Thus, diffusion along field lines is the mo st efficient transport, any ta ngling of field 
lines inside the gyroradius reduces the diffusion ( Colafrancesco fc Blasil 19981 ) . The electron 
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drift is hampered by pitch angle scattering with plasma waves, in particular with Alfven 
waves (IWentzell l Il974[ l. These waves are generated by fast electrons and protons that move 



with super-Alfven speed. Thus, collective effects in the plasma suppress the drift parallel 
to the field lines. Even if the precise treatment of electron diffusion is impossible since the 
actual state of the downstream plasma is not very well known, the collective effects seem to 
restrict the diffusion of suprathermal electrons to much smaller scales than those that we 
discuss in the context of radio relics, see Sec. El This allows us to assume that relativistic 
electrons purely advect with the plasma. Similar conclusions have been drawn in related 
studies that show, for example, that radio halos cannot be explained by the d iffusion of 
electrons from a central radio galaxy ( Jaffd 1977 : Dennisoiil IQSOl : Fuiita fc Sarazini.20oih . 



We now compute the total emission from the downstream region. The total emission per 
shock surface area A is given by the integral over the whole downstream region, the width 
of which is here parameterized by y 

"^dA d?^ = / /" dEn£;(i?,t)P(7L,a;obs). (29) 

In any given volume element, the electron spectrum depends on the time passed since the 
acceleration at the shock front due to cooling. If we assume that the suprathermal electrons 
advect with the plasma at constant speed, the distance to the shock front and the time 
passed since acceleration are related hy y = v^t. Hence, we can write 

~ I y ^cool 

6cooi et = Ccooi 7l — = 7^ — ry. 
Using the following substitution 

the integral Eq. ()29p becomes 

d^PKbs) ^^^^ r drr-(l-rr/)^-^Ffl^ ^C^^Ms). (31) 



dAdP ""7o Jr W 

The lower limit, Tmin, is in principle defined by the energy where the energy spectrum goes 
from thermal to suprathermal. However, as discussed above, the bulk of the radiation comes 
from the region with r ~ 0.03, which is by orders of magnitude above the value of Tmin. 
Thus, we are free to use as lower integration limit. The rapid decrease of F(l/r^ — » oo) 
ensures that the integrand is well-behaved. As a result, the integral depends only on the 
spectral index of the electron spectrum, s. 

Substituting the normalisation of the electron spectrum, Eq. (j2ip . the emission of a ho- 
mogenous volume Eq. (jTTj) . and Eq. (jnUj) to the total emission Eq. (|HT|). we find 



^ 2 ^sync ^ 1 ^t'^A H H Cy^Ud 

L^all — ITIqC —— X TT-eC/gpec 2 ^ 't' ~ ''^e ^/gpec ^sync L/,- — p:^ , 
O7- TJIqC ^cool ^cool 



where we have also used Vd = C^^/ud, see Eq. (fT^ . Now, it is useful to combine all factors 
that depend on the strength of the shock in Eq. ()31|) in 

(4.85 X 10"^)* (g - 1) /(7 - 1) qr^ - 1 Jvp 



2.07 X 10-9 q y gr7 (r - 1) 1,^,, ' 
where the normalisation constant dj, has been determined by the condition \E'(7V1 — * oo) = 1. 
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We note that there is also a weak dependency from the downstream plasma temperature via 
the integral /spec- However, for weak shocks, ^ 3, the expression \l/ decreases rapidly, see 
Fig.m while it is virtually constant with \E' ~ 1 for strong shocks. This implies that rather 
strong shocks with Mach numbers ^ 3 are necessary to produce any observable radio 
emission. 

After subsuming the terms governed by the strength of the shock under \E'(A^), we rewrite 
the total emission Eq. ^T\i 



du ^^^'^ \flG J \ Z/obs / C'cool C** 

= 6 4 X 10^4^ A_ ( '^"bs \ . . 

sHz Mpc2 10-4cm-3 0.05 Vl.4GHzy ^ ^ 

V7keV; /bcmbV . ( b'^ 

Note that the emission follows a power-law oc v^^^'^ and thus, for s = 2 it goes as oc 
Since ^!{M) ^ 1 and fii+"/V(5^MB + B^) ^ I for s = 2, the factor in front, 6.4 x 
]^q34 gj-gg-i Hz~^ Mpc~2, poses an upper limit for the synchrotron emission with respect to 
the Mach number of the shock and the strength of the magnetic field. 

If a radio relic is seen edge-on, its extension perpendicular to the shock front can be 
determined. From our model, we can compute at which distance, Uc-, the synchrotron emission 
in the downstream region decreases to a fraction (which we somewhat arbitrarily set to 30%) 
of its value immediately behind the shock front. We use such a rather high fraction since up 
to 10 % of the emission comes from electrons with r < 0.01, see Fig.El By the help of Eq. 
we can derive the corresponding Lorentz factor. For instance, in a region with a magnetic 
field of 4 /iG and with an observing frequency of 1.4 GHz, 10% of the synchrotron emission 
is originated by electrons with 7l < 100. The spectral electron density in this energy range, 
however, would be significantly diminished by Coulomb cooling, which is not included in 
our analysis. Thus, the distance where the emission diminishes to 30 % of the initial value 
seems to be a realistic estimate for the extension of the observable region. 

To obtain the width of the emission region, we derive first the dimensionless distance rjc, 
where rj has been introduced in Eq. ()30|) . The emission depends on the distance to the shock 
front, i. e. ?], via the integral J ... dr in Eq. ()3H) . We compute at which r/c this integral drops 
to 30 % of its value for r/ = 0. Since the integration depends on the slope of the electron 
spectrum, s, also rj^ does. We find a linear dependence, namely ?7c(s) = 4.8 x s — 1.5. Using 
the definition of rj we derive the effective width of the emission region 

*==cv;7* = ™""= {s^)\US' ^^^^ Tool^* 

This equation shows the dependence of the width of the relic, Uc-, on the magnetic field. In the 
regime of small magnetic field strengths, B <^ -Bcmb = 3.2/iG (for z = 0), i/c is proportional 
to 5^/2. Higher magnetic fields produce thicker relics because a higher magnetic field lowers 
the relevant Lorentz factor for synchrotron emission. In contrast, in a regime with strong 
magnetic field, B ^ -Bcmb? the Uc shrinks with increasing B-field according to yc oc B^^/"^ 
due to the synchrotron cooling itself. The width, Uc-, is maximal for i? ~ 2 /iG, see Fig.|Sl 
Also the strength of the shock affects the size of the synchrotron emitting area. Rather 
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weak shocks, ^ 3, are more extended, since the spectrum is steeper. For a spectral index 
of s = 3, electrons with smaller r are mainly responsible for the emission than for s = 2, 
see Fig. 01 Hence it takes longer to diminish the density of the relevant electron population. 
For small r, i. e., weak shocks and strong magnetic fields, the electron population is strongly 
reduced by Coulomb cooling, which is not included in our model. For those shocks the 
extension of the emission area is smaller than predicted here. 

We can also derive the time after which the emission of the downstream plasma has 
dropped to 30 % of its initial value from Eq. (jHHI) . 

tc — ~ n 

*-^coo1 

It has a maximum at S ~ 2 /iG, which amounts to ~ 1 Gyr for A4 = 4. 

The code for computing the synchrotron emission and the width of the emission region 
will be provided upon request. 



3 CLUSTER RADIO RELIC EXAMPLES 



3.1 AbelI-3667 

A textbook example of a merging cluster with strong shocks traveling into the cluster 
periphery is Abell3667 . This cluster is unique since it shows prominent radio emission 
( Rottgering et al.lll997l l which is presumab ly caused by a n outgoing shock wave. In X-rays, 
there is a double peak in t he cluster centre fehiba ta et al ]ll999) and a cold front indicating a 
fast moving substructure (jMarkevitch et al.. 2002 ). Both, indicates that the cluster suffered a 
merger recently. Even if the X-ray flu x from the cluster perip hery is not sufficient to identify 



a shock front at the position of relics, iRoettiger et al.l ()1999l ) showed that X-ray morphology 



and position of the relics is consistent with a merger of a 20 % substructure 1 Gyr ago. The 
spectral index of the north-west relic varies clearly: It steepens from a ~ 0.6 at the outer 
rim to 1.5 at its faint rim towards the cluster centre. This suggests that the relic is basically 
seen edge-on: The outer rim is populated by recently accelerated electrons, while towards 
the centre, the electrons population has aged. 

A similar steepening of the electron spectrum is expected in our model. As the plasma 
moves downstream and the electrons age, the radio spectrum steepens. However, the overall 
spectrum is a power-law, see Eq. even if locally the spectrum may steepen towards 
higher frequencies. Note that the power-law is only obtained if all contributions to the total 
emission are present: from the newly accelerated electron population at the shock front up 
to the old electrons that are too cold to emit significant synchrotron radiation. In order to 
apply our model to this relic, the downstream region should be homogeneous, i. e. it should 
not show adiabatic contraction i n this region, etc. . The n orth-west relic in Abell 3667 is quite 
close to such an ideal scenario. Rottgering et all (jl997l ) found an overall index of a = 1.1. 
They argued that there is no spectral steeping, what may indicate that we see plasma from 
the injection to the point where the emission has faded. Using Eq. ()32j) . we conclude that the 
initial slope in the electron spectrum is s = 2.2, which is consistent with a ~ 0.6, close to 
the shock surface (For a homogenous distribution of suprathermal electron with power-law 
spectral slope, s, the radio emission has the spectral index a = (s — l)/2). By combining 
Eq. (jSI), ()11|) and (fT^ . we can also infer the Mach number of the shock. We find that the 
Mach number of the shock related to the north-west relic in A 3667 is A^ = 4.7. 

We can use our expression for th e total emission to estimate the field strength in the 
downstream region. iRottgering et al.l (|l997|) derived the temperature and electron density 
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in the rehc region, 10® K and 10"'* cm~^, respectively. SarazinI ( 1999( ) argued that about 
1 % of the thermal energy in the ICM is found in relativistic electrons, with 7l > 300. 
Since our definition of includes also supra-thermal electron with 7l < 300, we choose a 
somewhat higher value of = 0.05. In order to compute th e total emission, we also need 
the area of the shock surface. From iRottgering et al.l (1997), Fig. 4, we estimate that the 
bulk of radio emission is extended over ~ 2.0 Mpc. Using these ICM parameters, we derive 
the radio emission as function of the shock strength and the downstream magnetic field, see 
Fig.ini The observed total emission is Pi. 4 = 4.1 x 10^^ ergs'^Hz"^. We have derived the 
Mach number of the shock by the spectral index of the radio emission. By comparing the 
observed radio power to the computed emission, see Fig. El we conclude that the strength of 
the magnetic field in the relic region is i? ~ 0.2 fiG. 

Since this relic is basically seen edge-on, its extension perpendicular to the shock front, 
i. e. radial with respect to the cluster centre, can be used to infer the magnetic field strength. 
In Sec. 12. 51 we have introduced a width i/c, which is the shortest distance between the shock 
front and the downstream position where the radio emission has decreased to 30% of its 
value at the shock front. The temperature of the downstream plasma and the Mach number 
of the shock determine the dimensionless quantity 77 and the downstream velocity, v^- For 
the parameters of A3667, we obtain 9.0 and 720kms~^, respectively. With a magnetic field 
of 0.2 fiG we derive by Eq. fjS^ that the width i/c amounts 0.3 Mpc, which is consistent with 
the observed extension of the relic. 

If the gas in the downstream region is substantially compressed, the relativistic electrons 
are re-energized by adiabatic compression. For a polytropic gas with 7 = 5/3, the gas density 
in the periphery of cluster scales with pgas oc (Infl + r/rJ/fr/rJ)'^/^, where r is the d istance 
from the cluster centre and is the scaling radius of the cluster (jAscasibar et al. 2006). 
Most relics are situated on the periphery of galaxy clusters. The radio relics in A3667, for 
example, are located at distances of ~ 2.6 Mpc on either side from the cluster centre. Given 
that the relics have a width of around 300 kpc and assuming that the downstream gas 
density follows the density profile of the cluster, the density increase along the rehc is no 
more than y 15%. Such a small cornpressi on has a moderate effect on the average electron 
energy, c.f. EnBlin fc Gopal-Krishnal (|2001 ). Flux conservation leads to i? oc p"^^^, hence the 



magnetic field can be enhanced maximally a few ten percent. Hence, we neglect the effects 
of adiabatic compression in the following discussion. 



3.2 AbelI-115 

Another example of a very extended radio relic has been found in Abell 155. Its size tan- 
gential to the cluster X-ray emission is more than 2 Mpc. In contrast, it is very thin in 
radial direction, a few hundred kpc at most. Compared to A3667, it has a much lower ra- 
dial/tangential extension ratio. Obviously, such a thin radio structure can only be achieved 
when the line-of-sight is tangential to the shock surface. Another prerequisite is that the ra- 
dio emission fades quickly while the plasma moves downstream. Considering the discussion 
above, this implies a low magnetic field strength in the downstream plasma. 

Similar to the procedure used for the north-west relic in A3667, we can infer the magn etic 
field strength in the emission region. We adopt the values given by IShibata et al. (|l999f ) for 



A115. The cluster tem perature in t he periphery is about 5keV, and the electron density of 
the order of 10~^cm~^. iGovoni et al.. (.2001.) investigated the radio emission and found a total 
emission of P1.4 = 1.9 x 10^^ ergs~^ Hz~^ with a slope of a = 1.1. We use those values and 
keep the other parameters to the values given above, e.g. the energy fraction in suprathermal 
electrons ^e- Since the spectral slope is the same as in A3667 the Mach number of the shock 
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is again = 4.7. The lower temperature in the downstream region also lowers the radio 
emission. However, the large difference to A3667 in radio power can only be achieved if the 
magnetic field is lower. We find that i? ~ 0.1 /iG reproduces the observed emission well. 

In the direction perpendicular to the shock front, the emission has dropped to 30% of 
its initial value at a distance of 0.1 Mpc from the shock front. This value is also consistent 
with the observations. 



3.3 AbelI-2256 

Abell2256 shows a radio halo and relic ( Clarke fc Enssliiil 2006f l. The relic is seen close to 
face-on, hence the spectral steepening is only mild. T he average slope is a = — 1-2, which 
indicates a shock with Mach number M. = 3.3. From Rephaeli fc Gruber ( 20031 ) we adopt 
a cluster temperature of 7.8 keV. The observed extension of the relic is about 1 Mpc, hence 



10 



31 



ergs 



Hz- 



is 



we assume a shock surface area of 1 Mpc^. The power, P1.4 = 3.6 x 
reproduced provided the magnetic field has a strength of -B ~ 0.3 fiG. 

This va lue for the magnet i c field is almost a magnitude lower as that in the central halo 



region (see IClarke fc Ensslh] (j200fi[) ) and indicates that the magnetic field decreases with 



increasing radius from the centre. 

If the radio relic is being produced by a giant shock wave, what has produced the radio 
halo in Abell 2256? Can the two phenomena be linked? With the magnetic field in the plasma 
also the relevant Lorentz factor for the radio e mission changes. If the m agnetic field in the 
halo region is about a few /iC as indicated by Clarke fc EnssliiJ ( 20061 ). the Lorentz factor 
of those electrons that produce the radio emission is 7L,haif ~ 300. In contrast, in the relic 
region the lower magnetic field implies that the Lorentz factor of the radio-emitting electrons 
is 7L,haif ~ 1000. At the field strengths in both relic and halo, we can assume that the cooling 
of the electron spectrum is mainly causes by inverse- Compton cooling. Since the cooling time 
is inversely proportional to 7l, the cluster core would emit synchrotron radiation three times 
longer. Hence it is conceivable that the halo and relic are caused by the same shock front. 



4 RADIO EMISSION VERSUS CLUSTER TEMPERATURE 

The example clusters discussed above indicate that shocks which produce presently observ- 
able radio emission have Mach numbers of about ~ 3 — 5. It is interesting to note that 
significantly stronger shocks show only moderately more radio emission if all other down- 
stream parameters are the same, as is shown in Fig.HJ This is due to the fact that strong 
shocks reach an asymptotic slope for the ele ctron spectrum. Much weaker shocks would have 



virtually no radio emission. iGabici fc Blasit (,2003.) found that it is quite unlikely to achieve 



shock strength of > 4 in a cluster merger. The combination of the paucity of shocks with 
Ai > 4 and the strong suppression of radio emission in weaker shocks may explain to some 
extent why only so few radio relics have been found. 

If we consider the magnetic field dependency in Eq. (j32|) . one striking feature is that 
for strong field strengths, B S> -Bcmb, the emission is constant. Hence, even an arbitrarily 
strong shock with an arbitrarily strong magnetic field would not exceed an emission of 
~ 6 X lO^^ergs"^ Hz~^, with standard parameters for temperature and density. Thus, there 
is a distinct upper limit for the radio emission. Intriguingly, our few examples support 
the assumption that the magnetic field in the ICM is correlated with its temperature. We 
found that A3667 and A2256 which have temperatures of about 8keV lead to magnetic field 
strengths of 0.2-0.3 fiG. In contrast, A115 with a temperature of 5keV shows a field strength 
of 0.1 /iG in the relic region. However, we restrict ourselves here to the custom assumption 



Radio signature from structure formation 15 



that a fraction, ^b, of the energy dissipated at the shock front is converted into magnetic 
fields, 



?2 



— = mpUeU. (34) 

On average, the magnetic field energy fraction is {C,b) = 6 x 10"'^ for our example clusters, 
see Tab-H Thus, only a very small fraction of the dissipated energy has to be converted. 
For A115 and A2256, the values of C,b niay be underestimated to some extent, since we have 
assumed that, at the location of the relics, the plasma has the average cluster temperature. 
Only for A3667, a temperature estimate for the relic region is available. However, we use 
Eq. (jH^ to compute the radio power as a function of the cluster temperature, see Fig.[7l 
solid line. Here we have assumed that the typical strength of a strong shock in the ICM is 
A4 = 4 and that the length scale of the shock fronts is 1 Mpc. 

Our estimate for the radio emission is close to the emission observed for radio halos. This 
may indicate that a significant part of the radio halo emission may be the result of an earlier 
merger shock. The observed radio power-temperature relation seems to be steeper than our 
result. This may be attributed to scaling relations that we have ignored in our analysis. 
It seems plausible that the linear size of radio halos depends on the cluster temperature, 
see Fig.lS] In average one may expect that the extension of halos is a fraction of the virial 
radius, r^ir which scales with the cluste r temperature by r^i^ oc M^/^ oc T^/^ (for cluster 
scaling relations see e. g. Eke et al. W9^. Including the linear size-temperature relation into 



our estimate for the radio power, we obtain a somewhat steeper dependency, see Fig.|7| The 
slope of the radio power is now dP/du oc 2"(i2+s)/4^ (j32p . where we 

have used 5^ oc T and A oc T. Hence, the strong temperature dependence is caused by 
the combination of several facts: With increasing cluster temperature, the electron density in 
the relevant Lorentz factor regime increases, the magnetic field scales with the temperature, 
and the average linear size of the emission increases with cluster temperature, i. e. its mass. 



5 COSMIC ACCRETION SHOCKS 



Numeri cal simulations su ggest that structure formation shocks may be subdivided into two 
classes f Rvu et al. 2003( ). Merger (or 'internal') shocks are caused by the supersonic mo- 
tion of infalling substru ctures, as seen, e.g., in A 3667 or in the 'bullet cluster' IE 0657-56 
" Markevitch et al.ll2nn2h . Typically, the Mach number of these shocks is rather low, A4 S 3, 



I^^rttaHWi^— ■/ ^ ±' ' ~ ^ 7 - ~ ~ ~ " ~ ~ ~ - -~ ---- - - - - _ y 

20031 ). since the cluster temperature reflects already the depth of the gravi- 



tational potential which also determines the maximal velocities in the course of the merger. 
When the shocks travel outwards into the cluster surroundings, the Mach number increases 
since the gas there has virtually never been shock-heated and is therefore much colder. Nu- 
merica l simulations indica t e that those accretion (or 'external') shocks reach Mach n umbers 
> 10 (|Miniati et al.ll2000l : IR vu et allEioi IPfrommer etHI 120061 : iHoeft et al.ll2006h . Basi- 
cally, clusters are envel oped by a sphere of accreti on shocks at a dist ance, rsh, of a f ew times 
the virial radius, r^i,, (iNagai fc Kravtsovi r2003: Rvu et al.l E)03: Ho eft et al.ll2006h . Due to 
the accretion shock, the temperature of the gas rises from ~ 10^ — 10^ K to almost the cluster 
temperature. To obtain a rough estimate for the temperature in the cluster periphery, we 



utilize the profile given bv lLoken et al. (|2002h 



T = 1.3 (T) 1 + 1.5 



-1.6 
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where (T) is the emission-weighted temperature of the cluster. With ~ 3rvir the temper- 
ature in the downstream region of accretion shocks is ~ 0.09 (T). Hence, we can determine 
the total emission caused by accretion shocks surrounding clusters of galaxies. For a cluster 
with (T) = 20 keV and rvir = 2 Mpc, we estimate the temperature in the downstream re- 
gion of the accretion shocks, namely 1.7 keV. Assuming that the electron density at such a 
large distance from the cluster centre amounts to 10~^ cm~^, we derive the magnetic field 
by Eq. we obtain 0.01 /iG. With the help of Eq. (jH^ . we compute the total emission of 
a sphere with 12 Mpc diameter. We obtain for a strong shock, Ai ^ 10, 

dP(1.4 GHz) 



du 



~ 2.2 X 10^^ erg s"^ Hz" 

-sh 



In reality, the accretion shock geometry is much more complex and projection effects have 
to be considered. Therefore, a significant fraction of the emission may come from a small 
region. Note that for such weak magnetic fields the emission region has only a small extension 
perpendicular to the shock front (~ 50 kpc). 

Upcoming radio telescopes such as Lofar and SKA will allow to explore the radio sky 
with signi ficantly impro v ed sen sitivity and resolution in the frequency range from 10 MHz 
to 1 GHz. Keshet et al. ( 2004[ ) inferred that the radio emission from structure formation 



shocks may become detectable with these telescopes. Hence, in near future we will have 
much stronger constraints for modeling the relativistic electron populations in clusters of 
galaxies. 

If we assume that 3 % of the emission comes from a rather compact region, we expect 
at 30 MHz a spot-like emission of ~ 4 x 10^^ ergs"^ Hz~^ (note P oc z/~*/^). If this cluster 
is located at 2; = 0.1 the fiux would amount to 1 mJy. This is close to the sensitivity of 
Lofar, hence accretion shocks with optimal conditions for radio emission will presumably 
be observed with sensitive radio telescopes in near future. 



6 SUMMARY 

Galaxy clusters grow either by mergers or by steady accretion, both of which produce large 
shock fronts. In particular, the extended radio relics {radio gischt) in the periphery of clus- 
ters are believed to trace these fronts. In analogy to supernova remnants, diffuse shock 
acceleration may produce a population of relativistic electrons that cool in the downstream 
region of the shock by inverse Compton and synchrotron losses. Here, we compute the total 
emission in the downstream region per unit area of the shock surface. The resulting analytic 
expression shows several interesting features: (i) the total spectrum is a power-law even if 
locally the spectrum may steepen towards higher frequencies. The overall slope is close to 
unity and depends solely on the slope of the initial electron spectrum which, in turn, is di- 
rectly related to the strength of the shock. Thus, the overall slope of the radio emission may 
allow us to infer the Mach number of the shock, (ii) the radio emission depends almost like 
a step-function on the shock strength. For Mach numbers < 3, only very little radiation is 
emitted while for > 10 the emission saturates, (iii) in the regime of weak magnetic fields, 
B <^ Bqmb, the total emission is proportional to B^+^Z^ while for strong magnetic fields the 
total emission is virtually independent of B. (iv) The width of the emission region, i. e. the 
extension perpendicular to the shock front, has a maximum extension for B ~ -Bcmb- 

We analysed the relics in A 115, A 2256, and A 3667 and found that the shocks have Mach 
numbers in the range of 3 to 5. Weaker shocks would lead to practically no radio emission. 
Furthermore, we can infer the magnetic fields in the emission regions. We find that for three 
sample relics the magnetic field is in the range from 0.1 to 0.3/iG. The extension of the relics 
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perpendicular to the shock front serves as an independent indicator for the field strength. 
The energy fraction stored in the magnetic field seems to be rather low, namely ~ 10~^. 

Assuming that the magnetic field energy density is on average proportional to the ther- 
mal one, we derive a radio emission-temperature relation. It agrees strikingly well with that 
found for radio halos. Moreover, if in the centres of galaxy clusters, i. e. where the radio 
halos reside, the magnetic fields is about a few /xG, the extension of the emission region 
perpendicular to the shock front would be several hundred kpc. Therefore, it is conceivable 
that at least part of the halo emission stems from shock waves which swept over the cluster 
centres within the last 1 Gyr. Using the value obtained for ^b, we have finally estimated the 
emission from cosmic accretion shocks. We inferred that at 30 MHz a radio power of 10^^ - 
10^°ergs~^ Hz"^ may be expected from regions where the accretion shock surface is parallel 
to the line-of-sight. A detailed computation of the cosmic radio background will be subject 
of a forthcoming study. 
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Table 1. Observed and derived parameters for our three model relics. For all computations we have assumed an electron 

density of no = 10"^ cm" ^ LLS indicates the longest linear extension of a relic. To determine the radio power of the relics we 
use Hq = 50kms— ^ Mpc— ^ in order to be consistent with computations in literature. 
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Figure 1. Entropy ratio, q = Sd/Su, compression ratio, r = pd/pu, and spectral index, s, of tlie electron spectrum generated 
by diffuse shock acceleration as a function of the Mach number. 
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Mach number 



Figure 2. Contours of the fraction of rclativistic electrons witli 7l > 300. Contour levels are at 10"**, 10-'', IQ-*, etc.. The 
energy fraction of suprathermal electrons is Jelec = 0.05. The initial cutoff in the spectrum is at e = lO'^. 
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Figure 3. Cumulative synchrotron emission. Computed according to Eq. 1271 assuming a power-law for the electron energy 
spectrum, nsir) oc r~*. 
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Figure 5. The distance at which the radio emission decreases to 30 % of the value at the shock front. 
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Figure 6. Radio power as a function of Mach number for different magnetic field strengths. We have marked the position of 
the radio power and the derived Mach number for Abell 3667. 
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T [keV] 

Figure 7. Radio power for an observing frequency of 1.4 GHz versus cluster temperature. The solid line indicates the radio 
power for a shock front with an area of 1 Mpc^, while the dashed line is computed for a temperature dependent ar ea, A ocT, see 
FigJ^ Radio and X-ray data for the clusters are taken from iGiovannini et al.. 1.1999.) , .Govoni et al. l20oI1) , and iGovoni et alJ 
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Figure 8. Longest linear size versus cluster temperature for radio halos shown in Fig.|7| The solid line indicates LLS = 
1 Mpc{T/7keV)i''2 rp 



